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Abstract 

A thermodynamically consistent density functional-perturbation theory is 
used to study the isostructural solid-to-solid transition which takes place in 
the hard sphere /attractive Yukawa system when the Yukawa tail is sufficiently 
short-ranged. A comparison with results for the square well potential allows 
us to study the effect of the attractive potential form on the solid-solid transi- 
tion. Reasonable agreement with simulations is found for the main transition 
properties as well as for the phase diagram evolution with the the range of 

the attractive potential. 
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I. INTRODUCTION 



The phase diagrams of colloidal suspensions have become a subject of growing interest 
in the last few years. Recent experimental results for mixtures of spherical colloidal par- 
ticles and nonadsorbing polymers show very interesting behaviour when the size of 
the polymer is reduced: The liquid-vapour critical temperature decreases and, eventually, 
it becomes smaller than the triple point temperature giving rise to a phase diagram where 
only the fluid and solid phases are stable and no liquid-vapour transition appears. A the- 
oretical understanding of this behaviour has been obtained from simplified models which 
substitute the polymer by an effective interaction which is attractive over distances larger 
than the colloidal particle diameter a. The origin of this attraction is an entropic effect || . 
However, the important point, as far as the aim of this paper is concerned, is that the range 
of this attraction is related to the polymer size and, therefore, reducing the polymer size 
means reducing the range of the attractive interaction between the colloidal particles. This 
mechanism gives the possibility of producing real "simple" systems with varying range of 
the interaction potential. Simple expressions (like the square-well potential or hard spheres 
with an attractive Yukawa tail) have been used for the resulting effective potential between 
colloidal particles in theoretical [HQ and computer simulation ||[7| studies of the liquid 
phase stability as a function of the range of the attractive interaction. Reasonable agree- 
ment between theory and simulation have been found (see figures [1-3] in ref. [§J) and the 
experimental results are qualitatively well understood. In summary, it is concluded that 
when the range of the attractive tail is smaller than « a/3 the liquid phase is not present in 
the phase diagram. In this case, when the gas density is increased at fixed temperature, the 
system reduces its free energy transforming itself into a solid before a van der Waals loop 
is generated. The sublimation line is then always above the liquid-gas coexistence curve in 
the density-temperature phase diagram and the system only shows solid and fluid as stable 
phases. 

Even more interesting phenomena take place if the range 5 of the attractive tail is further 
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reduced below a/3. Although direct experimental evidence does not exist so far, both 
simulations [[!]] and theory ]5|,|g|-|T0|j indicate that when 5 is small enough an isostructural solid- 



to-solid transition takes place between an expanded solid and a condensed one. This novel 
transition has some common features with the gas condensation, although its mechanism 
is different [||]. It is a first order transition that takes place between two phases with the 
same structure. It ends at a critical point as the temperature is increased and its location 
(in the density-temperature plane, for example) strongly depends on the range 5 of the 
attractive tail. However, some qualitative differences with respect to the condensation are 
rapidly observed. For example, the solid-solid transition is strongly asymmetric around 
its critical density and, contrary to the gas condensation, the qualitative way in which its 
critical temperature behaves with 5 crucially depends on the form of the attractive tail. In 
this way, T c remains practically constant as 5 is varied for the square well (SW) system and 
decreases with S for the hard sphere/attractive Yukawa potential (HSAY). 

In this paper we report theoretical predictions for the solid-solid transition of the HSAY 
potential, 

)oo , r < a 

(1) 
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obtained from our recently proposed density functional-perturbation theory ||11|| . In this 
work we will focus on extremely short-ranged attractive potentials (na > 25). As already 
pointed out, this theory gives reasonable agreement with simulations for smaller k values 
(na < 6) where the phase diagram shows its usual aspect with gas, liquid and solid stable 
phases. It is also able to predict reasonably the transition to a phase diagram without 
stable liquid phase for larger na values (see ref [|J). As shown below, our results for the 
solid-solid transition are also in qualitative agreement with simulations. Therefore, our 
density functional-perturbation theory gives a complete description of the phase diagram 
evolution with the range of the attractive interaction. 

The remainder of the paper is arranged as follows. In section II we briefly summarize our 
theory. We present and discuss our results in section III where we also include a comparison 



with simulations and with the predictions from other theoretical models |8|-|T0 



II. DENSITY FUNCTIONAL-PERTURBATION THEORY FOR SOLIDS 

The density functional formalism has been successfully applied to describe the solid 
phase of a hard sphere (HS) system in such a way that we have today very accurate density 
functional approximations for different properties of this system like, for example, its melting 
transition or its structure when it is placed in front of a hard wall, probably two of the 



strongest test for any theory of non- uniform classical systems | 12| . All these theories for 
non-uniform HS are based on some kind of mapping into a uniform system at some effective 
density. Different recipes to obtain the effective density from the physical one, p(r), give rise 
to different density functional approximations. However, this strategy will not work when 
it is applied to systems with more realistic potentials, like (p]), which include attractive 
interactions between particles. The reason is that the thermodynamic functions of the 
uniform fluid are not well defined inside the density gap corresponding to condensation 



13| , [T^| . It could be argued 0,0] that this is not an important problem when we are interested 



in cases where, like the ones we study in this paper, the liquid-gas transition is not present 
in the phase diagram because the range of the attractive tail is very small. However, there 
are reasons to believe that, even in these cases, the mapping strategy is problematic [15]. In 



effect, we would like to have a theory able to describe the complete phase diagram evolution 
with the potential range. It is obvious that this is not possible with the mapping strategy. 
Moreover, although the liquid-gas transition does not appear in the phase diagram when 
the potential range is small because it is preempted by the fluid-solid transition, the fluid 
free energy branch with the usual van der Waals loop corresponding to condensation does 
exist at low temperatures ||. Therefore, the mapping strategy will fail in this temperature 
region. 

An alternative is obtained from perturbation theories (PT) which have been successfully 



used in uniform fluids ]16[ . The standard perturbation scheme provides a density functional 
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approximation for the Helmholtz free energy, F [p(r)], of non-uniform systems. Its general 
expression is: 

F [p(r)] = F ref [p(r)] + U p [p(v)] . (2) 

F re f is the Helmholtz free energy of the reference system and U p is the perturbation contri- 
bution to the free energy. Up to first order, U p is given by 

U P [P( r )] = g / drdr 'pref( r > T ')M\ r - r 'l) ; ( 3 ) 

where <Pp( r ) is the perturbation and p re f(r, r ) is the pair distribution function of the reference 
system. This function is usually written in terms of g re /(r,r'), the extension of the radial 
distribution function (rdf) to non-uniform systems, as 

Prefix r ') = P( r )p( r ')9ref(r, r). (4) 

The reference system is now mapped into an equivalent HS system of diameter dns 
whose free energy is obtained from any of the accurate free energy functionals for non- 
uniform HS available in the literature. We use Tarazona model [17] in this work since it 



gives good results for the solid equation of state up to very high densities (~ lAd 3 HS ). This 
is a required condition to describe the solid-solid transition because the condensed solid 
density reaches values near the close packing limit when the temperature is low and/or the 
range of the attractive potential is very small. We do not need to specify a criterion to 
obtain the equivalent HS diameter dns from the reference potential in this case because the 
natural way to split the potential (]l|) into a reference and a perturbation parts already gives 
a hard sphere system of diameter a as reference. However, we need some approximation for 
the reference rdf g re f(r, r ) because we know very little of it in the non- uniform case. This 
problem is usually addressed making again some mapping into the uniform limit at some 
effective density p: 

g ref (r,r) w g u (\r - r'\,p), (5) 



where g u (\r — r'|, p) is the rdf of the uniform reference (HS in this case) system at density p. 
The effective density may in general depend on the position. It remains to give some recipe 
to obtain p. It is obvious that any recipe based on the local density p(r) will not work in 
the solid phase because p(r) reaches values which are many times larger than the maximum 
allowed density in a uniform fluid. A mapping like ([5]) has been used in PT for the Lennard- 
Jones against a wall flT8|J19H or for its melting transition |2(| using recipes for p which are 



the same or similar than those used for the free energy mapping. However, there are no 
reasons to think that the solid correlation structure and its free energy can be mapped into 
the uniform limit at equal or similar effective densities. It is clear that the solid free energy 
mapping is obtained using effective densities approximately equal to the mean solid density 
(this is the result from current density functionals for HS). However, the effective density we 
need in the rdf mapping (at least as far as the perturbation energy calculation is concerned) 
is expected to be very small (and then g u should be almost a step function). The reason 
is that most of the correlation structure of p^ in (H) is already implicitly included in the 
density product of equation (|4]). To take into account this idea we have recently proposed 
an alternative based on imposing an exact sum rule to obtain p (see ref. Jll|]). This is the 
local compressibility equation: 

/*>Wr,r)-l] = -l + ^. «0 

where p is the chemical potential. Precisely, p is defined through equation (|) substituting 
the non-uniform rdf g re /(r, r') by its uniform limit g re f(\ r ~ r '| ; p(r)) and solving for p(r). A 
simpler version of this theory may be obtained using a constant effective density p which is 
obtained from the exact global compressibility equation [TJJ] : 



J drp(r) J drp(r')\g ref (r,r) - 1] = -1 + k B Tpx T , (7) 

where xt is the solid isothermal compressibility and p its mean density. Using this recipe 
we obtain p < 0.08. This is consistent with our previous intuition on g u and also means that 
the mean field approximation obtained using a step function for g u should be a reasonable 



approximation for the perturbation contribution to the free energy of the solid. Notice also 
that Eqn. (|6|) or (|7p reduces to the well known compressibility equation in the uniform 
limit. Therefore, our theory reduces to the standard perturbation theory in this limit and, 
in this way, we have a theory which treats both phases on the same footing. This is crucial 
since we are then able to determine consistently the relative stability of one phase with 
respect to the others as it is required to study the corresponding phase transitions and, 
therefore, to obtain the phase diagram. Our theory is also able to continuously go from 
one phase to the other if the local compressibility equation (||) with a non-constant effective 
density is used. This allows a consistent study of the solid-liquid interface. Moreover, it does 
not suffer from the same problems than the non-perturbative approximations (as discussed 
above) and, therefore, it is able to describe the phase diagram for arbitrary values of the 
range of the attractive potential. Then we can study the complete phase diagram evolution 
using always the same theory. This is an important point in order to give predictions for 
the relevant values of ko where the phase diagram suffers interesting qualitative changes 
like, for example, the disappearance of the liquid phase or the appearance of the solid-solid 
transition. 



III. RESULTS AND DISCUSSION 

We have studied the phase diagram of the HSAY potential ([!]) for ko = 25, 33, 40, 
50 and 67 using the PT of the previous section in its simplest form, i.e. using a constant 
effective density p obtained from the global compressibility equation (|7]). The election of 
these particular values of kg has been motivated by the existence of simulation results for 
those values 0. We restrict ourselves to study the face-centered-cubic (fee) structure for the 
solid phase. As usual, we write the solid local density, p(r), as a sum over the lattice sites of 
normalized gaussian peaks. The free energy functional (0)is minimized with respect to the 
gaussian width parameter. The standard common tangent method is then used to obtain the 
coexisting densities of the different phases. The same theoretical approach was used in our 
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previous work |4j] on the same potential where we studied the liquid phase extinction as kg 



is increased and the attractive potential becomes more and more short-ranged. Therefore, 
those results together with the ones we are presenting in this section show a complete and 
consistent picture of the phase diagram evolution. 

Our results are shown in figures 1-5 as continuous lines. The theoretical results for tran- 
sitions between metastable states are shown by dashed lines. We also show the simulation 
results of Bolhuis et al. (dotted lines). Note that the simulation results for the solid-liquid 
transition for the cases kg = 40, 50, and 67 have not been reported. The main characteris- 
tics of the solid-solid transition summarized in the Introduction are qualitatively reproduced 
by the theory. Contrary to the SW case [||,[7| , the critical point temperature T c depends now 
on the range of the attractive tail. In fact, both simulations and theory indicate that T c is 
almost a linear function of k. We can give arguments to understand this difference: The 
attractive energy per particle (in k^T units) may be reasonably approximated, when the 
range of the attractive potential is small enough, by 



where pr 1 (r) is the normalized gaussian density peak centred at the nearest neighbour shell 
distance Ri and n n is the number of nearest neighbours (12 in the fee lattice). A further 
simplification can be introduced considering the low values of p given by (|7|) (as discussed 
in the previous section) and, therefore, using the mean field approximation. Then g u is 



From this last equation it is easy to see that the form of (p p (r) plays an important role. 
For the SW at low solid densities (i.e. those corresponding to the expanded solid) the 
gaussian peak p Rl (r) is, roughly speaking, practically outside the attractive well. Then the 
attractive energy is very small. On the contrary, in the condensed solid side the gaussian 
peak is completely inside the well. Then, since the potential is constant and the gaussian 




(8) 



approximated by a step function and Eqn (§) can be written as 




(9) 
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peak is normalized, the attractive energy is ~ — 6e//c#T. The important point is that this 
result depends on the temperature only. It does not depend on the well width nor on the 
density. Therefore, the critical point temperature is not expected to depend on this two 
parameters. However, the same argument does not apply to the HSAY potential. The main 
reason is that the form of the potential does not allow to say when the density peak is inside 
the well and when it is outside. The angular integration in (§) can be straightforwardly 
done to obtain the attractive energy in terms of a single radial integral whose result will 
obviously depend on k. At a given density its absolute value decreases when k increases. 
The solid- solid critical temperature should then be lower when the range of the attractive 
potential is reduced. It is interesting to stress the dependence of the critical temperature 
behaviour with the well width on the potential form. This is qualitatively different from the 
liquid-vapour critical temperature case. This last temperature is always reduced when the 
well width is reduced, irrespective of the potential form. The reason resides in the different 
transition mechanisms. The solid-solid transition takes place when a particle is able to 
see the attractive wells due to its neighbours. The precise moment when this happens 
depends on the well form. A liquid particle, however, always sees the attractive wells of the 
surrounding particles. Its attractive energy is given (in the mean field approximation) by 
the total integral of the potential and this always depends on the potential width and form. 

We turn now to a more detailed comparison with the simulation results. The theoret- 
ically predicted solid-solid critical temperature shows a practically linear dependence with 
k with a slope in agreement with simulations. However, our result for T c is systematically 
above the simulations with an overestimation around 13.5%. Some authors have attributed 
this overestimation to the mean field nature of the theory, in the same way that mean field 
approximations overestimate the liquid- vapour critical temperature ||. On the contrary, we 
believe that the origin of this discrepancy between theory and simulation is different. As 
discussed in the previous section, the mean field approximation for the solid phase takes into 
account most of the solid correlation structure. Moreover, we have obtained an overestima- 
tion around 59% in our prediction for the solid-solid critical temperature of the SW potential 

Q 



H using exactly the same theory. It is difficult to believe that correlations together with 
the different form of the potential may be responsible for this enormous difference. Instead, 
we think that the solid-solid transition is the result of the very delicate balance between 
the different contributions to the free energy of the solid phase ||. In fact, the final result 
crucially depends on particular details of the model. A further confirmation is the fact that 
we obtain /c^T c /e ~ 2.7 for the SW potential using the mean field approximation (i. e. 
using a step function for g re f) for the attractive energy together with Tarazona model for 
the HS reference system. This result should be compared with fcgT c /e ~ 2.0 obtained by 
Likos et al. using the very same mean field approximation for the attractive energy but 



using Denton and Ashcroft functional for HS [21]. Moreover, Likos et al. have studied the 



solid-solid transition of the SW potential using a non-perturbative approximation for the 
full potential |10[ . While they get improved results for the critical temperature as compared 



with their previous mean field calculation || , they obtain now a critical temperature which 
depends on the potential range, in opposition to simulation and mean field results. 

Fig. 1-5 show that our theoretical result for the critical solid- solid density is slightly 
smaller than the simulation result. However, this underestimation reduces when k increases. 
We also overestimates the value of k where the solid-solid transition emerges in the phase 
diagram (see figures 1-3). It is interesting to note that the predictions from the same theory 
for these two properties in the SW case are practically identical to those of simulations ||. 
As already discussed, the form of the potential makes now the transition mechanism less 
decisive in order to fix the precise localization of the transition. There is not an absolute 
criterion to say when a gaussian density peak is inside the attractive wells due to its nearest 
neighbours and when it is outside. This means that, contrary to the SW case, the results 
for the relevant densities (like the critical and triple points densities) will now depend on 
the model details because the nearest neighbour distance is only function of the mean solid 
density. However, this problem becomes less important when k increases because the solid- 
solid transition moves to higher densities, near the close packing limit, where everything is 
geometrically fixed. 

in 



In conclusion, we find a solid-solid transition whose coexistence curve is, as far as its gen- 
eral form is concerned, in good agreement with simulation (see figures 3-5) but shifted in the 
density-temperature plane. The results for the relevant properties (critical temperature and 
density, etc.) depend on model details because the solid-solid transition is the consequence 
of a delicate balance between the hard core energy and the attractive energy. At the same 
time, this means that the form of the attractive potential plays an important role which 
determines some qualitative aspects like the dependence of the critical temperature with the 
range of the potential. We have shown that our perturbation theory, outlined in section II, 
is able to describe reasonably the phase diagram of the HSAY system for the range of the 
attractive potential ranging from "usual" values (where the phase diagram shows liquid-gas, 
liquid-solid and vapour-solid transitions) to extremely short-ranged potentials (where the 
liquid-vapour condensation is not present) and even shorter-ranged potentials where the 
novel isostructural solid-solid transition emerges in the phase diagram. To our knowledge, 
this is the only consistent theory able to perform it successfully. 



ACKNOWLEDGMENTS 

We thank Dr. C.N. Likos for sending us a preprint of ref. [[Hj prior to publication. This 
work was supported by the Direction General de Investigation Cientffica y Tecnica of Spain, 
under Grant PB91 - 0090. 



1 1 



REFERENCES 

[1] F. Leal Calderon, J. Bibette and J. Biais, Europhys. Lett. 23, 653 (1993). 

[2] P.N. Pusey, W. C. K. Poon, S. M. Ilett, and P. Bartlett, J. Phys.: Condens. Matter, 6, 
A29 (1994). 

[3] E.J. Meijer and D. Frenkel, Phys. Rev. Lett. 67, 1110 (1991). 

[4] L. Mederos, and G. Navascues, J. Chem. Phys., 101, 9841 (1994). 

[5] C. Rascon, G. Navascues, and L. Mederos Phys. Rev. B, 51, 14899 (1995). 

[6] E.J. Meijer and D. Frenkel, J. Chem. Phys., 100, 6873 (1994). 

[7] P. Bolhuis, M. Hagen, and D. Frenkel, Phys. Rev. E, 50,4880 (1994). 

[8] C.F. Tejero, A. Daanoum, H.N.W. Lekkerkerker, and M. Baus, Phys. Rev. Lett. 73, 
752 (1994). 

[9] C.N. Likos, Zs.T. Nemeth, and H. L6wen, J. Phys. Condens. Matter, 6, 10965 (1994). 
[10] C.N. Likos, and G. Senatore J. Phys. Condens. Matter, in press. 

[11] L. Mederos, G. Navascues, P. Tarazona, and E. Chacon, Phys. Rev. E. 47, 4284 (1993). 
See also L. Mederos, G. Navascues, and P. Tarazona, Phys. Rev. E, 49, 2161 (1994). 

[12] See, for example, Fundamentals of Inhomogeneous Fluids, edited by D. Henderson 
(Dekker, New York 1992). 

[13] A. de Kuijper, W.L. Vos, J.-L. Barrat, J. -P. Hansen, and J. A. Schouten, J. Chem. Phys. 
93, 5187 (1990). 

[14] A. Kyrlidis and R.A. Brown, Phys. Rev. A. 44, 8141 (1991). 

[15] C. Rascon, G. Navascues, and L. Mederos J. Phys. Condens. Matter, in press. 

[16] J. P. Hansen and I.R. McDonald, Theory of Simple Liquids (Academic, NY, 1986). 

1 9 



[17] P. Tarazona, Phys. Rev. A 31, 2672 (1985). 

[18] Z. Tang, L.E. Scriven, and H.T. Davis, J. Chem. Phys. 95, 2659 (1991). 
[19] S. Sokolowski and J. Fischer, J. Chem. Phys. 96, 5441 (1992). 
[20] A. Kyrlidis and R.A. Brown, Phys. Rev. E. 47, 427 (1993). 
[21] A. Denton and N.W. Ashcroft, Phys. Rev. A 39, 4701 (1989). 



1 3 



FIGURES 

FIG. 1. Phase diagram of the hard sphere/attractive Yukawa potential in the tempera- 
ture-density plane for na = 25. The solid lines are the predictions of our theory. The dashed 
lines are also predictions from our theory but for transitions between metastable states. The 
dotted lines are the simulation results of Bolhuis et al. 



FIG. 2. As figure 1 for no = 33. 



FIG. 3. As figure 1 for kg = 40. 



FIG. 4. As figure 1 for no = 50. 



FIG. 5. As figure 1 for no = 67. 
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